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On using visibility correlations to probe the HI 
distribution from the dark ages to the present epoch I: 
Formalism and the expected signal 
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ABSTRACT 

Redshifted 21 cm radiation originating from the cosmological distribution of neutral 
hydrogen (HI) appears as a background radiation in low frequency radio observations. 
The angular and frequency domain fluctuations in this radiation carry information 
about cosmological structure formation. We propose that correlations between visi- 
bilities measured at different baselines and frequencies in radio-interferometric obser- 
vations be used to quantify the statistical properties of these fluctuations. This has 
an inherent advantage over other statistical estimators in that it deals directly with 
the visibilities which are the primary quantities measured in radio-interferometric ob- 
servations. Also, the visibility correlation has a very simple relation with the power 
spectrum. We present estimates of the expected signal for nearly the entire post- 
recombination era, from the dark ages to the present epoch. The epoch of reionization, 
where the HI has a patchy distribution, has a distinct signature where the signal is 
determined by the size of the discrete ionized regions. The signal at other epochs, 
where the HI follows the dark matter, is determined largely by the power spectrum 
of dark matter fluctuations. The signal is strongest for baselines where the antenna 
separations are within a few hundred times the wavelength of observation, and an op- 
timal strategy would preferentially sample these baselines. In the frequency domain, 
for most baselines the visibilities at two different frequencies are uncorrelated beyond 
IS.v ~ 1 MHz, a signature which in principle would allow the HI signal to be easily 
distinguished from the continuum sources of contamination. 

Key words: cosmology: theory - cosmology: large scale structure of universe - diffuse 
radiation 



1 INTRODUCTION 

One of the major problems in modern cosmology is to de- 
termine the distribution of matter on large scales in the 
universe and understand the large scale structure (LSS) 
formation . The last decade has witnessed phenomenal 
progress, particularly on the observational front where 
CMBR anisotropics (eg. Spergel et al. 2003) and galaxy red- 
shift surveys (eg. Tegmark et al. 2003a, 2003b) in conjunc- 
tion with other observations have determined the parame- 
ters of the background cosmological model and the power 
spectrum of density fluctuations to a high level of preci- 
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sion. These observations are all found to be consistent with 
a ACDM cosmological model with a scale invariant primor- 
dial power spectrum, though there are hints of what could 
prove to be very interesting features beyond the simplest 
model. 

Another very interesting problem, on smaller scales in 
cosmology, is the formation of the first gravitationally bound 
objects and their subsequent merger to produce galaxies. 
Observations of quasar spectra show the diffuse gas in the 
universe to be completely ionized at redshifts z < 5 (Fan et 
al. 2002). Understanding the reionization process and its re- 
lation with the first generation of luminous objects produced 
by gravitational collapse is one of the challenges facing mod- 
ern cosmology. 

Observations of the redshifted 21 cm HI radiation pro- 
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Figure 1. This figure schematically shows the evolution of HI 
(shaded in the figure) through the entire history of the universe. 
Hydrogen decouples from the CMBR at z ~ 1000 when electrons 
and protons combine, for the first time, to produce HI. The sub- 
sequent era, until the formation of the first luminous objects at 
z ~ 30, is referred to as the dark ages or the Pre-reionization era. 
The HI, in this era, is nearly uniformly distributed with small 
fluctuations which trace the dark matter. The first luminous ob- 
jects reionize the HI. The reionization is initially confined to small 
bubbles surrounding the luminous objects. The bubbles of ionized 
HII gas grow until they finally overlap and the whole universe is 
reionized by z ~ 6. The HI distribution in the Reionization era 
or the Gray Age is patchy. Only the HI in dense collapsed ob- 
jects survive the reionization. The surviving HI clouds which are 
seen as DLAs in quasar absorption spectra are shown as points 
in the figure. The large scale distribution of these clouds is as- 
sumed to follow the dark matter. Also shown is the CMBR which 
propagates through the intervening HI to the observer at z = 0. 



vide an unique opportunity for probing the universe from 
the dark ages when the hydrogen gas first decoupled from 
the CMBR to the present epoch. Such observations will al- 
low us to probe in detail how the universe was reionized. 
In addition, redshifted HI observations will allow the large 
scale distribution of matter in the universe to be studies over 
a very wide range of redshifts. The possibility of observing 
21 cm emission from the cosmological structure formation 
was first recognized by Sunyaev and Zeldovich (1972). Later 
studies by Hogan & Rees (1979) and Scott & Rees (1990) 
consider both emission and absorption against the CMBR. 

The cosmological evolution of the HI gas is shown 
schematically in Figure 1. The 21 cm emission has been 
perceived as a very important tool for studying the epoch 
of reionization, and the expected signature, in 21 cm, of 
the heating of the HI gas and its subsequent reionization 
has been studied extensively (Madau, Meiksin & Rees 1997 
Gnedin & Ostriker 1997; Shaver et al. 1999; Tozzi et al. 2000 
Iliev et al. 2002; Iliev et al. 2003; Ciardi & Madau 2003 
Furlanetto,Sokasian. & Hernquist 2003; Gnedin & Shaver 
2003; Miralda-Escude 2003; Chen & Miralda-Escude 2004). 
The timing, duration and character of the events which led 
to the reionization of the universe contains an enormous 
amount of information about the first cosmic structures and 
also has important implication for the later generations of 
baryonic objects. This process is not very well understood 
(eg. Barkana & Loeb 2001), but there are several relevant 
observational constraints. 

The observation of quasars at redshift z ~ 6 which show 
strong HI absorption (Becker et al. 2001) indicates that at 
least 1% of the total hydrogen mass at z ~ 6 is neutral (Fan 
et al. 2002), and the neutral mass fraction decreases rapidly 



at lower redshifts. This is a strong indication that the epoch 
of reionization ended at z ~ 6 . 

Observations of the CMBR polarization, generated 
through Thomson scattering of CMBR photons by free elec- 
trons along the line of sight, indicates that the reionization 
began at a redshift z > 14. On the other hand, the observed 
anisotropics of the CMBR indicate that the total optical 
depth of the Thomson scattering is not extremely high, sug- 
gesting that reionization could not have started at redshift 
much higher than about 30 ((Spergel et al. 2003). 

A third constraint comes from determinations of the 
IGM temperature from observations of the Lycv forest in 
the z range 2 to 4 which indicates a complex reionization 
history with there possibly being an order unity change in 
the neutral hydrogen fraction at z < 10 (Theuns et al. 2002; 
Hui& Haiman 2003). 

The possibility of observing the HI emission at z < 6, 
the post-reionization era, has also been discussed extensively 
particularly in the context of the Giant Meterwave Radio 
Telescope (GMRT) (Subramanian & Padmanabhan 1993; 
Kumar, Padmanabhan & Subramanian 1985; Bagla, Nath 
and Padmanabhan 1997; Bharadwaj, Nath & Sethi 2001; 
Bagla & White 2002). These observations will probe the 
large scale structure formation at the redshift where the HI 
emission originated. 

Loeb & Zaldarriaga (2003) have recently proposed that 
observations of the angular fluctuations in the HI absorp- 
tions against the CMBR, originating in the z range 30 to 
200 corresponding to the dark age of the universe, would 
allow the primordial power spectrum to be probed at an 
unprecedented level of accuracy. Bharadwaj & Ali (2004) 
have investigated this in detail including the effect of red- 
shift space distortion and gas temperature fluctuations. 

In this paper we discuss the nature and the magnitude 
of the HI signal expected across the entire redshift range 
starting from the dark age of the universe at z ~ 1000 
through the epoch of reionization to the present era. Our dis- 
cussion is in the context of observations using radio interfer- 
ometric arrays to study the redshifted HI signal. Such obser- 
vations are sensitive to only the angular fluctuations in the 
HI signal. The quantity measured in radio interferometric 
observations is the complex visibility. Typically, the visibil- 
ities are used to construct an image which is then analysed. 
An image of the HI distribution, while interesting, is not the 
optimal tool for quantifying and analysing the redshifted HI 
radiation. The quantities of interest are the statistical prop- 
erties of the fluctuations in the redshifted HI radiation, and 
an image, if at all useful, is only an intermediate step in 
the process of extracting these statistical quantities. Here 
we propose that correlations between the complex visibili- 
ties measured at different baselines and frequencies be used 
to directly quantify the angular and frequency domain clus- 
tering pattern in the HI radiation, completely doing away 
with the need for making an image. 

In Section 2 of this paper we develop the formalism for 
calculating visibility correlations, bringing out the relation 
between this quantity and various properties of the HI dis- 
tribution. The fluctuations, with angle and frequency, of the 
redshifted HI radiation can be directly related to fluctua- 
tions in the properties of the HI distribution in space. The 
mapping is a little complicated due to the presence of pe- 
culiar velocities which moves around spectral features along 
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the frequency direction. This effect, not included in most 
earlier calculations, is very similar to the redshift space dis- 
tortion in galaxy redshift surveys (Kaiser 1987) and has been 
included here (also Bharadwaj, Nath & Sethi 2001; Bharad- 
waj & Ali 2004). 

The possibility of using visibility correlations to study 
the HI distribution has been discussed earlier in the re- 
stricted context of HI emission from z < 6 by Bharadwaj 
& Sethi (2001), Bharadwaj & Pandey (2003) and Bharad- 
waj & Srikant (2004) who have developed the formalism and 
used it to make detailed predictions and simulations of the 
signal expected at the GMRT. 

Morales & Hewitt (2003) (MH hereafter) have discussed 
the use of visibility correlations in the context of detecting 
the HI signal from the epoch of reionization using LOFAR. 
They also address the issue of extracting the signal from the 
foreground contamination. 

Zaldarriaga, Furlanetto & Hernquist (2003) (ZFL here- 
after) have proposed the use of the angular power spectrum 
to quantify the fluctuations in the redshifted HI radiation, in 
direct analogy with the analysis of the CMBR anisotropies, 
As noted by them, though this analogy is partially true in 
that both the redshifted HI radiation and the CMBR are 
background radiations present at all directions and frequen- 
cies, there is a major difference in that observations at differ- 
ent frequencies probe the HI at different distances whereas 
we expect exactly the same CMBR anisotropies at all fre- 
quencies. The multi-frequency angular power-spectrum in- 
troduced by ZFL allows the cross-correlations in the angular 
fluctuations in the HI signal at different frequencies to be 
quantified. 

The cross correlation between the HI signal at different 
frequencies is expected to decay rapidly as the frequency 
separation is increased, and as noted by several earlier au- 
thors (eg. Shaver et al. 1999; DiMatteo. et al. 2002; Gnedin 
& Shaver 2003; DiMatteo. et al. 2004) this holds the possibil- 
ity of allowing us to distinguish it from the contaminations 
which are expected to have a continuum spectrum and be 
correlated across frequencies. ZFL discuss this issue focusing 
on the properties of the contaminants which could possibly 
limit the ability to extract the HI signal, and show that 
it should be possible to extract the signal provided point 
sources with flux > O.lmJy can be identified and removed. 

In this paper we build up on the earlier work in this 
field. In particular, we generalize the formalism developed in 
Bharadwaj & Sethi (2001) so that it can be used for both, 
HI emission and absorption against the CMBR, thereby ex- 
tending its scope to the entire redshift range starting from 
the epoch when hydrogen first recombines to the present. 
The formalism reported here presents progress over the ear- 
lier work (MH and ZFL) on at least two counts. First, it 
incorporates the effects of redshift-space distortions, ignored 
in all earlier works except those by Bharadwaj and collab- 
orators (mentioned earlier). As noted by Bharadwaj & Ali 
(2004), this is a very important effect and not including it 
changes the signal by 50% or more. Second, our formalism 
shows, explicitly, how the the correlations in the HI signal 
decay with increasing frequency separation Av. The Av de- 
pendence, though implicit in the formalism of MH and ZFL, 
has not been explicitly quantified by these authors. A third 
point, which in our opinion is also a major advantage of our 
formalism over the angular power spectrum (ZFL), is that it 



deals directly with visibilities which are the primary quan- 
tity measured in radio- interferometric observations. Further, 
the visibility correlations have a very simple relation to the 
three-dimensional power spectrum of HI fluctuations, the 
quantity of prime interest in HI observations. It may be 
noted that though the visibility-correlations and angular 
power-spectrum are equivalent descriptors of the redshifted 
HI observations, the former is simpler to calculate and inter- 
pret, it being related to the HI power spectrum through an 
exponential function instead of spherical Bessel functions. 

We next present a brief outline of the paper. In Section 
3, we trace the evolution of the cosmological HI from the 
dark ages to the present era, focusing on epochs which are 
of interest for 21 cm observations. In Section 4 we present the 
results for the expected HI signal from the different epochs, 
and in Section 5 we summarize and discuss our results. 

Finally, we bring to the notice of readers two recent 
papers (Morales 2004) and (Santos, Cooray & Knox 2004), 
submitted while this paper was being revised, addressing 
respectively the power spectrum sensitivity and the Av de- 
pendence of of the epoch of reionization HI signal. 



2 THE FORMALISM FOR HI VISIBILITY 
CORRELATIONS 

The story of neutral hydrogen (HI) starts at a redshift 
z ~ 1000 when, for the first time, the primeval plasma 
has cooled sufficiently for protons and electrons to com- 
bine. The subsequent evolution of the HI is very interest- 
ing, we shall come back to this once we have the formu- 
las required to calculate the HI emission and absorption 
against the CMBR. The propagation of the CMBR, from the 
last scattering surface where it decouples from the primeval 
plasma to the observer at present, through the intervening 
neutral or partially ionized hydrogen is shown schematically 
in Figure 1. Along any line of sight, the CMBR interacts 
with the HI through the spin flip hyperfine transition at 
f e = c/A e = 1420 MHz (A e = 21 cm) in the rest frame of the 
hydrogen. This changes the brightness temperature of the 
CMBR to 



T(r) = T 7 e- T + T s (l- e - T ) 



(1) 



where r is the optical depth, T 7 is the background CMBR 
temperature and T s is the spin temperature defined through 
the level population ratio 

- = ^exp(-7VT s ). (2) 

no go 

Here no and ni refer respectively to the number density 
of hydrogen atoms in the ground and excited states of the 
hyperfine transition, go = 1 and g\ = 3 are the degeneracies 
of these levels and T» = 0.068 A" = h p v e /kB, where hp is 
the Planck constant and fcs the Boltzmann constant. 

The absorption/emission features introduced by the in- 
tervening HI are redshifted to a frequency v = f e /(l + z) 
for an observer at present. The expansion of the universe, 
and the HI peculiar velocity both contribute to the redshift. 
Incorporating these effects, the optical depth at a redshift z 
along a line of sight n is (Appendix A) 
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pHI 

pk 



(1 + z) dv 
H(z) dr v 



(3) 



where r„ and H(z) are respectively the comoving distance 
and the Hubble parameter, as functions of z or equivalently 
v, calculated using the background cosmological model with 
no peculiar velocities. Also, Phi/ph is the ratio of the neutral 
hydrogen to the mean hydrogen density, and v is the line of 
sight component of the peculiar velocity of the HI, both 
quantities being evaluated at a comoving position x = r„n. 
The optical depth (eq. 3) is typically much less than unity. 
The quantity of interest for redshifted 21 cm observa- 
tions, the excess brightness temperature redshifted to the 
observer at present is 



ST b (n,z) = 



T(t) - T-, (T 3 - T 7 )t 



(4) 



1+2 1 + Z 

It is convenient to express the excess brightness tem- 
perature as STb(n, z) — T(z) x r]ni{n,z) where 



T(z) =4.0mK(l + z) 



» b h 2 
0.02 



0.7\ H„ 



V h ) H 



(z) 



(5) 



depends only on z and the background cosmological param- 
eters, and 



Tjm(n,z) - 



Phi 

PH 



1- 



T 

_ 



1- 



(1 + z) dv 
H(z) dr 



(6) 



is the "21cm radiation efficiency" (1 — T 7 /T s ) phi/ph intro- 
duced by Madau, Meiksin & Rees (1997), with the extra 
velocity term arising here on account of the HI peculiar ve- 
locities. We refer to i/hi as the "21 cm radiation efficiency 
in redshift space". This varies with position (x = rvn) and 
redshift, and it incorporates the details of the HI evolution 
and the effects of the growth of large scale structures. We 
also introduce ?7Hi(k, z), the Fourier transform of r/m(x, z), 
defined through 



??Hi(n,2) = 



d A k 



— i k ■ r u n - 



fjm(K z) 



(7) 



(2tt) 3 

and the associated three dimensional power spectrum 
Pm(k, z) defined through 

(J7Hi(k, z) 77 HI (k', z)) = (2^) 3 S%(k - k') Pki(k, z) (8) 

where 8% is the three dimensional Dirac delta function. 

When discussing radio interferometric observations, it 
is convenient to use the specific intensity I v instead of the 
brightness temperature. In the Raleigh- Jeans limit SI V = 
[2k B /\l(l + z) 2 ] x ST b = I v x -qm where 



I v = 2.5 x 10^ 



Jy fn b h 2 \ /0.7\ H 



0.02 



(x) 



H(z) 



(9) 



The quantities T and I v are determined by cosmological 
parameters whose values are reasonably well established. In 
the rest of the paper, when required, we shall use equations 
(5) and (9)_with {Q m o, fixo, ^bh 2 , h) = (0.3,0.7,0.02,0.7) to 
calculate T and I v . The crux of the issue of observing HI 
is in the evolution r/fn which is largely unknown. We shall 
discuss a possible scenario in the next section, before that we 
discuss the role of 77m in radio interferometric observations. 
We consider radio interferometric observations using an 
array of low frequency radio antennas distributed on a plane. 
The antennas all point in the same direction m which we 



take to be vertically up wards. The beam pattern A{9) quan- 
tifies how the individual antenna, pointing up wards, re- 
sponds to signals from different directions in the sky. This 
is assumed to be a Gaussian A(9) = e~ e ^ e ° with 8q <gC 1 
i. e. the beam width of the antennas is small, and the part of 
the sky which contributes to the signal can be well approx- 
imated by a plane. In this approximation the unit vector n 
can be represented by n = m + 9, where 9 is a two dimen- 
sional vector in the plane of the sky. Using this SI V can be 
expressed as 

j3 1 



Sh(n)=h 



d A k 

0^)3 



-ir u (fci|+kj_.0)- 



rjm(k,z) 



(10) 



where k» — k- m and k^ are respectively the components of 
k parallel and perpendicular to m. The component kj_ lies 
in the plane of the sky. 

The quantity measured in interferometric observations 
is the complex visibility V(\J, v) which is recorded for ev- 
ery independent pair of antennas at every frequency chan- 
nel in the band of observations. For any pair of antennas, 
U = d/A quantifies the separation d in units of the wave- 
length A, we refer to this dimensionless quantity U as a 
baseline. A typical radio interferometric array simultane- 
ously measures visibilities at a large number of baselines and 
frequency channels. Ideally, each visibility records a single 
mode of the Fourier transform of the specific intensity distri- 
bution I v {9) on the sky. In reality, it is the Fourier transform 
of A(9) I v {9) which is recorded 



V{\J,v) = 



d 2 9A{9) I v {6) e - l27rU ' e . 



(11) 



The measured visibilities are Fourier inverted to de- 
termine the specific intensity distribution I v (9), which is 
what we call an image. Usually the visibility at zero spacing 
U = is not used, and an uniform specific intensity distribu- 
tion make no contribution to the visibilities. The visibilities 
record only the angular fluctuations in I v {ff). If the spe- 
cific intensity distribution on the sky were decomposed into 
Fourier modes, ideally the visibility at a baseline U would 
record the complex amplitude of only a single mode whose 
period is 1/U in angle on the sky and is oriented in the direc- 
tion parallel to U. In reality, the response to Fourier modes 
on the sky is smeared a little because of A(9), the antenna 
beam pattern, which appears in equation (11). 

We now calculate the visibilities arising from angular 
fluctuations in the specific intensity excess (or decrement) 
SI V (9) produced by redshifted HI emission (or absorption) 
against the CMBR. Using eq. (10) in eq. (11) gives us 



V(U, v) = /„ 



^3 o(U - ^k ± )r; H i(k, z)e-^^ (12) 



where a(U) the Fourier transform of the antenna beam 
pattern A(ff) 



A(9) = 



,2 Tr -i2itU-6 tTT \ 

due a\U) . 



(13) 



For a Gaussian beam A{9) — e~ 6 ' e o , the Fourier transform 
also is a Gaussian a(U) = ty9q exp [— tt 2 9qU 2 ] which we use 
in the rest of this paper. 

In equation (12), the contribution to V(\J, v) from dif- 
ferent modes is peaked around the values of k for which 
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kx = 27rU/r„. This is so because the visibility V{\J,v) re- 
sponds to angular fluctuations of period 1/U on the sky, 
which corresponds to a spatial period r v /U at the distance 
where the HI is located. It then follows that 77m (k, 2) will 
contribute to a visibility U(U, v) only if the component of 
k projected on the plane of the sky satisfies kx — 2tt\J jr v . 
The effect of the antenna beam pattern is to introduce a 
range in kx of width | Akx |~ 2/r v 6o to which each visibil- 
ity responds. 

We use eq. (12) to calculate (V(U, v)V*(U' , v + Ai/)), 
the correlation expected between the visibilities measured 
at two different baselines U and U , at two frequencies v 
and v + Av which are slightly different ie. Av jv <C 1. The 
change v — » v + Av will cause a very small change in all the 
terms in eq. (12) except for the phase which we can write as 
e -ik||(rv+TvAi,) wllere r ' v _ dr v /du. Using this the visibility 
correlation is 

<V(U, v) V\V , v + Av)) =P v J^L a (V- gk x ) 



(U 



^kx)fti(k)e lfc ii r - A " 
Ait 



(14) 



The visibilities at U and U will be correlated only if there 
is a significant overlap between the terms a(U — f^-kx) 

and a*(U — Tp-kx) which are peaked around different val- 
ues of kx- It then follows that visibilities at two different 
baselines are correlated only if | U — U |< l/(n9o) and 
the visibilities at widely separated baselines are uncorre- 
cted. To understand the nature of the visibility correla- 
tions it will suffice to restrict our analysis to U = U , 
bearing in mind that the correlation falls off very quickly 
if I U — U |> l/(7T#o)- A further simplification is possible 
if we assume that U 3> l/2n8o for which we can approx- 
imate the Gaussian a(U) with a Dirac delta function us- 
ing a 2 (\ U (r„/27T)k x |) ~ (2n 3 9 2 /rl)8 2 D (k ± - (2ir/r„)U) 
whereby the integrals over kx in equation (14) can be eval- 
uated. We then have 



(V(U,v)V*(U,v+Av)) 



4ri 



dk l} PhiQs.) e l 



'(15) 



where k = /ciim + (2-7r/r„)U. We expect Pm(k) to be an 
even function of fcii. It then follows that the imaginary part 
of the visibility correlation is zero , and the real part is 



(y{\J,v)V{\5,v + Av)) = 



■Ivi 



(16) 



dk\\ -PHi(k) cos(fc||r^A^) 



Considering the case Av = first, we see that the 
visibility correlation (V(U, v)V*(XJ, v)) is proportional to 
Pm(k) integrated over k\\ . This directly reflects the fact that 
a baseline U responds to all Fourier modes k whose projec- 
tion in the plane of the sky matches kx = (2n/r v )\J. A 
mode k = fcym + (2n/r v )\J satisfies this condition for ar- 
bitrary values of k\\ which is why we have to sum up the 
contribution from all values of fcy through the integral. Fur- 
ther interpretation is simplified if we assume that Pm(k) 

is an isotropic, power law of k — « k 2 + (2ttU /r v ) 2 ie. 



Pui(k) oc k n . The integral converges for n < — 1 for which wc 
have {V(TJ,v)V*(TJ, v)) oc U 1+n and it is a monotonically 



decreasing function of U, ie. the correlations will always be 
highest at the small baselines. This behaviour holds in much 
of the range of r v and U which we are interested in. In gen- 
eral Pni(k) will not be an isotropic function of k, but the 
qualitative features discussed above will still hold. 

The correlation between the visibilities falls as Av is in- 
creased. To estimate Av where the visibility correlation falls 
to zero, we assume that Pni(k) oc k n with n < 0. The inte- 
gral in eq. (16) will have a substantial value if the power 
spectrum falls significantly before the term cos(fci|r„ Av) 
can complete one oscillation, and the integral will vanish 
if the power spectrum does not fall fast enough. Using the 
fact that the power spectrum will change by order unity 
when APhi/Pui ~ nAk«/{2'nU/r v ) ~ 1 and that the cos 
term completes a full oscillation when Ak\\r v Av ~ 2ti, 
we have an estimate Av ~ r v /(r v U) for the value of 
Av beyond which the visibilities become uncorrelated. The 
point to note is that the fall in the visibility correlation 
(U(U, v)V*(\J, v + Av)) with increasing Av depends on the 
combination r v jr v which in turn depends only on 2 and the 
cosmological parameters. This holds the possibility of us- 
ing HI observations to probe the background cosmological 
model over a large range of 2. 

It may be noted that we have assumed infinite frequency 
resolution when developing the formalism. Actual observa- 
tions will have frequency channels of finite width. We briefly 
discuss the implications of this in Section 5. 

In the next section we shall consider a model for the 
evolution of the HI, and use this to make predictions for the 
visibility correlations from redshifted 21 cm radiation. 



3 THE EVOLUTION OF THE HI 

The evolution of the HI is shown schematically in Figure 
1. We identify three different eras namely Pre-reionization, 
Reionization and Post-reionization and discuss these sepa- 
rately. 



3.1 Pre-reionization 

We resume the story of HI at the redshift 2 ~ 1000 where it 
decouples from the CMBR and the universe enters the Dark 
Age. Subsequent to this the HI gas temperature T g is main- 
tained at T~f by the CMBR which pumps energy into the 
gas through the scattering of CMBR photons by the small 
fraction of surviving free electrons. This process becomes in- 
effective in coupling T g to T 7 at 2 ~ 200. In the absence of 
external heating at 2 < 200 the gas cools adiabatically with 
T g oc (I + 2) 2 while T 7 oc (I + 2). The spin temperature T s is 
strongly coupled to T g through the collisional spin flipping 
process until 2 ~ 70. The collisional process is weak at lower 
redshifts, and T s again approaches T 7 . This gives a window 
in redshift 30 < 2 < 200 where T 3 < T-, and rjm < 0. The HI 
in this redshift range will produce absorption features in the 
CMBR spectrum. Loeb & Zaldarriaga (2003) have proposed 
that observations of the angular fluctuations in this signal 
holds the promise of allowing the primordial power spectrum 
of density perturbations to be studied at an unprecedented 
level of accuracy. The nature of this signal was studied in 
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detail in Bharadwaj & Ali (2004), incorporating the effects 
of peculiar velocities and gas temperature fluctuations. 

We assume that all the hydrogen is neutral, and that 
the fluctuations in the HI density Apni trace the dark 
matter fluctuations A ie ApHi(k, z) = pHA(k,z). Th HI 
peculiar velocities are also determined by the dark mat- 
ter, and the Fourier transform of [(1 + z) / H (z)](dv / dr) 
in equation (6) is — /(£7 m )p 2 A(k, z), where p is the co- 
sine of the angle between k and the line of sight n, and 
/(fi m ) « 0.° m 6 + i [l - ifi m (l + On)] in a spatially flat 
universe (Lahav et al. 1991). The fluctuations Apm also 
cause fluctuations in the spin temperature AT S which we 
quantify using a dimensionless function s(z) defined such 
that AT S — sT s Apm/pm- Using these in equations (6) and 
(7) we have the Fourier transform of the "21 cm radiation 
efficiency in redshift space" 



Vm( 



(.-!)(. + /,'v 



A(k,*) 



(17) 



which we use to calculate Phi, the power spectrum of the 
radiation efficiency in terms of the dark matter power spec- 
trum P(k) 



Pm(k,«) = 



£)(l + //*') + £« 



T, 



T s 



P(k, 



(18) 



Observations of the angular fluctuations in the HI radia- 
tion using either the visibility correlations or other meth- 
ods like the angular power spectrum will determine Phi (eq. 
16). These observations will directly probe the matter power 
spectrum provided the evolution of (1 — T~ t /T s ) and s Tj/T s 
are known. This has been investigated in detail in Bharad- 
waj & Ali (2004) and we use the results here. The quantities 
(1— Ty/T,,) and s T y /T 3 also determine the frequencies where 
we expect a substantial signal. 

We use eq. (18) to calculate the visibility correlations 
expected from the Pre-reionization era. 



3.2 Reionization 

The Dark Age of the universe ended, possibly at a redshift 
z ~ 30, when the first luminous objects were formed. The 
radiation from these luminous objects and from the subse- 
quently formed luminous objects ionized the low density HI 
in the universe. Initially only small spherical regions (Strom- 
gren sphere) surrounding the luminous objects are filled 
with ionized HII gas, the rest of the universe being filled 
with HI (Figure 1). Gradually these spheres of ionized gas 
grow until they finally overlap, filling up the whole of space, 
and all the low density gas in the universe is ionized. Here 
we adopt a simple model for the reionization process, The 
model, though simple, is sufficiently general to reasonably 
describe a variety of scenarios by changing the values of the 
model parameters. 

We assume that the HI gas is heated well before it is 
reionized, and that the spin temperature is coupled to the 
gas temperature with T s ^> T 7 so that (1 — T 7 /T s ) — » 1. It 
then follows that r/ui > (eq. 6) ie. the HI will be seen in 
emission. Also, the HI emission depends only on the HI den- 
sity and peculiar velocity. Like in the Pre-reionization era, 
we assume that the hydrogen density and peculiar velocities 
follow the dark matter, the only difference being that a frac- 
tion of the volume fv is now ionized. We assume that non- 




Figure 2. This shows the evolution of the mean neutral frac- 
tion Shi with the frequency of the redshiftcd HI emission for the 
reionization model discussed in the text. 



overlapping spheres of comoving radius R are completely 
ionized, the centres of the spheres being clustered with a 
bias b c > 1 relative to underlying dark matter distribution. 
This model is similar to that used by Zaldarriaga, Furlan- 
etto & Hernquist (2003) in the context of HI emission, and 
Gruzinov & Hu (1998) and Knox at al. (1998) in the con- 
text of the effect of patchy reionization on the CMBR. There 
is a difference in that we have assumed the centers of the 
spheres to be clustered, whereas Zaldarriaga et al. (2003) 
assumed them to be randomly distributed. One would ex- 
pect the centres of the ionized spheres to clustered, given 
the fact that we identify them with the locations of the first 
luminous objects which are believed to have formed at the 
peaks of the density fluctuations. 

The fraction of volume ionized fv, the mean neutral 
fraction xhi and mean comoving number density of ionized 
spheres tihi are related as fv = 1 — Shi = (47rP 3 /3)nHi- 
Following Zaldarriaga et al. (2003), we assume that 

XHl{z) = — J- ry-7-r (19) 

1 +cxp((z — z )/Az) 

with zo = 10 and Az = 0.5 so that 50% os the hydrogen 
is reionized at a redshift z — 10 (Figure 2). The comoving 
size of the spheres R increases with z so that it is 3 h~ 1 Mpc 
at z — 10, and R increases in such a way so that the mean 
comoving number density of the ionized spheres is constant 
at nm = 4.4 x 10~ 3 ft 3 Mpc~ 3 . 

In this model the HI density is pm(x, z) — Ph(1 + 
8) [l — ^2 a ^(1 x — Xa /^)] wnere 6 is the dark matter fluc- 
tuation, a refers to the different ionized spheres with centers 
at x a , and 6(y) is the Heaviside step function defined such 
that 9{y) = 1 for < y < 1 and zero otherwise. We then 
have 



r/Hi(x, z) — 



1 + 6- 



1 + z dv 
H(z)dr 



£< 



R 



.(20) 



We have already discussed the Fourier transform of the first 
part of eq. (20) which is (27r) 3 <5 3 3 (k) + (1 + p 2 )A(k, z), and 
the Fourier transform of the second part is (27r) 3 5f ;) (k) — 
(47rP 3 /3)W(fcP)X; a e lk ' x ° where W(y) = (3/y 3 )[sin(y) - 
ycos(y)] is the spherical top hat window function. The term 
y] e lkx ° is the Fourier transform of the distribution of the 
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centers of the ionized spheres which we write as ^ e !kx a _ 
n H i[(27r) 3 c5o(k) + A P (k, z) +& c A(k,«)], where A P (k,z) is 
the fluctuation in the distribution of the centers arising due 
to the discrete nature of these points (Poisson fluctuations) 
and the term fc c A(k, z) arises due to the clustering of the 
centers. These two components of the fluctuation are inde- 
pendent and (Ap(k)Ap(k')) = (27r) 3 (5|,(k - k') x l/n ffl . 
Convolving the Fourier transform of the two terms in eq. 
(20), and dropping terms of order A 2 and AAp we have 

r?Hi(k,2) = [x m (l + fv 2 )-b c f v W{kR)] A(k,z) 

- f v W{kR)A P (k,z) (21) 

This gives the power spectrum of the "21 cm emission effi- 
ciency in redshift space" to be 



Pm(k,z) 



[x m (l + f/J. 2 
fvW 2 (kR) 



b c f v W(kR)] 2 P(k, 



z) + 



riHi 



(22) 



The first term which contains P(k) arises from the cluster- 
ing of the hydrogen and the clustering of the centers of the 
ionized spheres. The second term which has 1/nni arises due 
to the discrete nature of the ionized regions. 

Our model has a limitation that it cannot be used 
when a large fraction of the volume is ionized as the ion- 
ized spheres start to overlap and the HI density becomes 
negative in the overlapping regions. Calculating the fraction 
of the total volume where the HI density is negative, we 
find this to be fyjl under the assumption that the centers 
of the ionized spheres are randomly distributed. We use this 
to asses the range of validity of our model. We restrict the 
model to z > 10 where fv < 0.5, and the HI density is neg- 
ative in less than 12.5% of the total volume. The possibility 
of the spheres overlapping increases if they are highly clus- 
tered, and we restrict b c to b c = 1.5 throughout to keep this 
under control. 

We use equation (22) to calculate the visibility correla- 
tions during the Reionization era. 

3.3 Post-reionization 

All the low density hydrogen is ionized by a redshift z ~ 6 
and HI is to be found only in high density clouds (Figure 1), 
possibly protogalaxies, which survive reionization. Observa- 
tions of Lyman-Q absorption lines seen in quasar spectra 
have been used to determine the HI density in the redshift 
range 1 < z < 3.5. These observations currently indicate 
Qgas(z), the comoving density of neutral gas expressed as a 
fraction of the present critical density, to be nearly constant 
at a value Q, gas (z) ~ 10~ 3 for z > 1 (Peroux et al. 2001). 
The bulk of the neutral gas is in clouds which have HI col- 
umn densities greater than 2 x 10 20 atoms/cm 2 (Peroux et al. 
2001, Storrie-Lombardi, McMahon & Irwin 1996, Lanzetta, 
Wolfe & Turnshek 1998). These high column density clouds 
are responsible for the damped Lyman-Q absorption lines 
observed along lines of sight to quasars. The flux of HI 
emission from individual clouds (< 10/iJy) is too weak to 
be detected by existing radio telescopes unless the image of 
the cloud is significantly magnified by an intervening cluster 
gravitational lens (Saini, Bharadwaj & Sethi 2001). 

For the purposes of this paper, we assume that the HI 
clouds trace the dark matter. It may be noted that we do not 



really expect this to be true, and the HI clouds clouds will, 
in all probability, be biased with respect to the dark mat- 
ter. However, our assumption is justified given the fact that 
we currently have very little information about the actual 
distribution of the HI clouds, 

Converting £l ga s to the mean neutral fraction Shi = 
pHi/pH = £l g as/Qb gives us xm = 50£lgash 2 (0.02/£l b h 2 ) or 
xhi = 2.45 x 10~ 2 . We also assume T s ^> T 7 and hence we 
see the HI in emission. Using these we have 



P HI (k,z) = xl 1 (l + fu 2 ) P(k,z) 



(23) 



The fact that the neutral hydrogen is in discrete clouds 
makes a contribution which we do not include here. This 
effect originates from the fact that the HI emission line from 
individual clouds has a finite width, and the visibility corre- 
lation is enhanced when Av is smaller than the line-width 
of the emission from the individual clouds. Another impor- 
tant effect not included here is that the fluctuations become 
non-linear at low z. Both these effects have been studied 
using simulations (Bharadwaj & Srikant 2004), The simple 
analytic treatment adopted here suffices for the purposes of 
this paper where the main focus is to estimate the magni- 
tude and the nature of the HI signal over a very large range 
of redshifts starting from the Dark Ages to the present. 



4 RESULTS 

In this section we present predictions for the visibility cor- 
relation signal from the different eras. Equation (16) allows 
us to calculate the visibility correlation by evaluating just 
a single integral. This is valid only if U 3> (2n6o)~ , failing 
which we have to use eq. (14) which involves a three dimen- 
sional integral. Bharadwaj & Pandey (2003) have compared 
the results calculated using equations (14) and (16), and find 
that they agree quite well over the range of U of interest. 
Here we use eq. (16) throughout. 

In equation (16), it is necessary to specify #o which is 
the size of the beam of the individual antennas in the array. 
It may be noted that 9o ~ 0.6 x 6>fwhm- The value of 9q will 
depend on the physical dimensions of the antennas and the 
wavelength at which the observations are being carried out. 
We have assumed that #o = 1° at 325MHz and scaled this 
as 6o oc A to determine the value at other frequencies. This 
value for 9q is appropriate for the GMRT which can operate 
at frequencies down to 150MHz. The visibility correlations 
scale as # 2 , and it is straightforward to scale the results 
presented here to make visibility correlation predictions for 
other radio telescopes. 

The other input needed to calculate visibility correla- 
tions is PHi(k), the power spectrum of the "21 cm radiation 
efficiency in redshift space" . In the models for the three eras 
discussed in the previous section, the only quantity which is 
left to be specified is the dark matter power spectrum P(k) 
which we take to be that of the currently favoured ACDM 
model with the shape parameter Y — 0.2, normalized to 
as — 1 at present. We use equations (18), (22) and (23) 
for Pfn(k) in the pre-reionization, reionization and post- 
reionization eras respectively. 

The visibility correlation (V(U, v)V*(U, v + Av)) re- 
sponds to the power spectrum at Fourier modes k > kmin = 
(2n/r v )U. Figure 3 shows the part of the power spectrum 
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Figure 3. This shows the linear power spectrum of dark matter 
density fluctuations used to calculate the expected HI visibility 
correlations. The power spectrum is shown at the present epoch. 
A baseline U at a frequency v will probe the power spectrum at 
all Fourier modes k > k m i„ = (2n/r^)U. The values of k m i n arc 
shown at different frequencies for U = 100. This figure can be 
used to determine the length-scale probed by any baselines U. 
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Figure 5. This shows the correlation between the visibilities 
V(U, v) and V(U, v + Ai/) expected for the same baseline U at 
two slightly different frequencies. The result are shown for differ- 
ent values of U (shown in the figure) for the 32MHz band. 



rapidly, falling by an order of magnitude by U = 1000. Also, 
with varying frequency the signal falls off quite fast beyond 
~ 45 MHz. Frequencies smaller than 25 MHz will be severely 
affected by the ionosphere and have not been shown. The 
visibility correlations as a function U have a very similar 
behaviour at all the frequencies, with just the amplitude 
of the signal changing with frequency. The shape of these 
curves directly reflect the power spectrum at the relevant 
Fourier modes. 

The visibility correlation (V(U, v)V*(XJ, v+Av) for the 
same baseline and different frequency is shown for 32 MHz 
in Figure 5. The behaviour is similar at other frequencies. 
The correlation falls rapidly as Av is increased. The value 
of Av where the visibilities become uncorrelated scales as 
~ r v /{r' v U) and it is around 0.8 MHz for U = 100. 



Figure 4. This shows how the visibility correlation 
(V(U,u)V*(U,v)) varies with baseline U. The results arc 
shown for the different frequencies shown in the figure. 



which will be probed at different frequencies by the visibility 
correlations at the baseline U — 100. The point to note is 
that a typical radio interferometric array with a large range 
of baselines and frequencies will provide many independent 
measurements of the power spectrum at different epochs and 
in overlapping bins of Fourier modes. This will allow the evo- 
lution of the HI and the dark matter power spectrum to be 
simultaneously studied at a high level of precision. 



4.1 Pre-reionization 

The signal here will be absorption features in the CMBR 
spectrum. The visibility correlations record the angular fluc- 
tuations in this decrement. The results for Av = are shown 
in Figure 4. The signal is maximum at a frequency ~ 32 MHz 
where the perturbations in the HI density are most efficient 
in absorbing the CMBR flux (Bharadwaj & Ali 2004). The 
signal strength at any fixed frequency does not vary with 
U in the range of baselines 10 to 100, after which it falls 



4.2 Reionization 

The HI signal here is in emission. The results for the visibil- 
ity correlation (V(U, v)V*(XJ, v)) are shown as a function of 
U for different frequencies in Figure 6. At v < 110 MHz the 
hydrogen is largely neutral, and the visibility correlations 
trace the dark matter power spectrum. The ionization picks 
up at v = 120 MHz (Figure 2) where ionized bubbles fill 
up ~ 20% of the universe. The centers of these bubbles are 
clustered with a bias b c — 1.5 with respect to the underly- 
ing dark matter distribution. The presence of these clustered 
bubbles reduces the signal at small baselines U < 400. At 
130 MHz around 50% of the universe is occupied by ion- 
ized bubbles. Here the clustering of the bubbles or of the 
underlying dark matter is of no consequence, and the vis- 
ibility correlations are due to the discrete nature (Poisson 
distribution) of the ionized bubbles. The visibility correla- 
tion falls drastically for the large baselines which resolve out 
the bubbles, and the signal at U > 1000 is due to the cluster- 
ing of the HI which follows the dark matter. The model for 
the reionization breaks down when a large fraction of the 
universe is ionized, and we do not present results beyond 
130 MHz. 

An interesting feature when the visibility correlation 
is dominated by the contribution from individual bubbles 
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Figure 6. This shows how the visibility correlation 
(V (U , v)V * (U , v)) varies with baseline U. The results arc 
shown for the different frequencies shown in the figure. 





l.4e 


-07 


-4 


1.2e 


-07 
















\e 


-07 


A 












<l 


8e 


-08 














;j 


6e 


-08 








* 






> 


4e 


-08 








y 


?.e. 


-08 


> 






V 









REIONIZATION 

130 MHz 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Av MHz 



Figure 7. This shows the correlation between the visibilities 
V(U, v) and V(U,v + Af) expected for the same baseline U at 
two slightly different frequencies. The result are shown for differ- 
ent values of U (shown in the figure) for the f30MHz band. 



Figure 8. This shows how the visibility correlation 
(V(U, u)V*(\J, v)) varies with baseline U. The results are 
shown for the different frequencies shown in the figure. 



4.3 Post-reionization 

The predictions for {V(U, v)V*(\J, v)) are shown as a func- 
tion of U for a variety of frequencies in Figure 8. At all 
frequencies the correlations are nearly independent of U in 
the range 10 to 100, and then fall by nearly an order of 
magnitude by U — 1000. The HI in this era has been as- 
sumed to trace the dark matter, and the shape of the curve 
showing the visibility correlation as a function of U is de- 
cided by the dark matter power spectrum at the relevant 
Fourier modes. The correlations increase with frequency, a 
reflection of the fact that the matter power spectrum grows 
with time. It may be noted that the signal also depends on 
Qgas which we have assumed to be a constant over redshift, 
and the evolution of the signal probes both the growth of 
the power spectrum and the evolution of ilgas. At a fixed 
baseline, the correlation {V(\J, v)V*(\J, v + Av)) falls with 
increasing Av. The value of Av where the visibilities be- 
come uncorrelated scales as Av ~ r^/r^U . The predictions 
for v = 500 MHz are shown in Figure 9 where the visibility 
at U = 100 becomes uncorrelated at Av ~ 1.3 MHz. The 
behaviour is similar at the at the other frequencies. 



is that the fall in {V(\J,v)V*(\J,v + Av) with increasing 
Av is decided by the ratio Rjr v . This follows from the fact 
that the power spectrum Phi is now proportional to W (kR) 
which falls rapidly for k 3> 1/R, and the oscillating integral 
in eq. (16) cancels out when Av > R/r v . The value of Av 
where the visibilities become uncorrelated depends only on 
the size of the bubbles R and does not change with U. 

It should be noted that studies of the growth of ion- 
ized bubbles (eg. Furlanetto, Zaldarriaga & Hernquist 2004) 
show them to have a range of sizes at any given epoch dur- 
ing reionization. This will possibly weaken some of the ef- 
fects like the sharp drop in visibility correlations (Figure 6) 
and the Av dependence (Figure 7) calculated here assuming 
bubbles of a fixed size. 



5 SUMMARY AND DISCUSSION. 

We have proposed that correlations between the visibilities 
measured at different baselines and frequencies in radio- 
interferometric observations be used to quantify and inter- 
pret the redshifted 21 cm radiation from cosmological HI. 
This statistics has an inherent advantage over other quanti- 
ties which have been proposed for the same purpose in that 
it deals directly with visibilities which are the primary quan- 
tities measured in radio-interferometric observations. There 
is a further advantage in that the system noise in the vis- 
ibilities measured at different baselines and frequencies is 
uncorrelated. 

The visibility correlations, we show, have a very simple 
relation to the power spectrum (equation 16), the latter be- 
ing the statistical quantity most commonly used to quantify 
the cosmological matter distribution. Also, there is a di- 
rect relation between the dimensionless baseline U and the 
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Figure 9. This shows the correlation between the visibilities 
V(U, v) and V(U, v + Av) expected for the same baseline U at 
two slightly different frequencies. The result are shown for differ- 
ent values of U (shown in the figure) for the 500MHz band. 



Fourier mode probed by the visibilities at that particular 
baseline. 

We have made predictions for the signal expected over 
the entire evolution of the HI starting from the dark ages 
to the present epoch. The magnitude of the signal depends 
on do ~ 0.6 x Ofwhm the size of the beam of the individ- 
ual antenna elements in the array. Our predictions are for 
0o = 1° at 325MHz which is appropriate for the GMRT. At 
other frequencies we have used 80 oc A in making our pre- 
dictions. The results presented here can be easily used to 
make predictions for other telescopes using the fact that the 
visibility correlations scale as 8q. The shape of the curves 
showing (V(U, v)V*(\J, v + Av)) as functions of U or Av 
have no dependence on #0 or any other parameter of the 
individual antennas. Finally, a reminder that our analysis 
assumes 80 <C 1 where the curvature of the sky may be 
neglected. The analysis will be more complicated for an ar- 
ray of small antennas where the individual elements have a 
very wide beam, but possibly the qualitative features of this 
analysis will hold. 

One of the salient features of our analysis is that for all 
the eras the visibility correlation signal from redshifted HI 
is maximum at small baselines. In the situation where the 
signal arises mainly from the large scale clustering of the 
HI whose distribution follows the dark matter the signal is 
nearly constant up to U ~ 100 after which the signal falls 
rapidly dropping by an order of magnitude by U ~ 1000. In 
the situation when the process of reionization is underway 
and the HI distribution is patchy, the behaviour of the visi- 
bility correlation is decided by the size of the ionized regions. 
In our model for reionization where there are spherical ion- 
ized bubbles of comoving radius R — 3 ft _1 Mpc, the visibil- 
ity correlation is a constant to U ~ 400, the larger baselines 
are not sensitive to the ionized spheres. Details of the reion- 
ization model aside, the baseline dependence of the visibility 
correlations is sensitive to size of the ionized patches and it 
should be possible to extract this information. Also, an op- 
timal observational strategy would preferentially sample the 
small baselines where the signal is largest. 



We next shift our attention to the correlation between 
the visibilities at the same baseline but at two different fre- 
quencies with separation Av. We see that the correlations 
fall off very fast as Av is increased. In the situation where 
the visibility correlation signal is from the large scale grav- 
itational clustering of the HI the value of Av where the 
visibilities become uncorrelated scales as r v /(r v U) and it is 
around 1 MHz at U = 100 at all frequencies. In the situa- 
tion where the signal is from ionized patches, the value of 
Av where the visibilities become uncorrelated is decided by 
R/r„ and it is around 0.4 MHz independent of the value of 
U . This feature should allow us to determine from the obser- 
vations whether it is the large scale gravitational clustering 
of the HI or the ionized patches which dominates the visi- 
bility correlation. The latter is expected to dominate during 
the epoch when reionization is underway. It may be noted 
that the value 0.4 MHz quoted here is specific to the model 
of reionization which we have adopted. Though the value 
is model dependent, we expect that the qualitative features 
discussed here should continue to hold even in a more gen- 
eral situation. 

The main challenge in observing cosmological HI is to 
extract the HI signal from the various contaminants which 
are expected to swamp this signal. The contaminants include 
Galactic synchrotron emission, free-free emission from ion- 
izing halos (Oh & Mack 2003), faint radio loud quasars (Di- 
Matteo. et al. 2002) and synchrotron emission from low red- 
shift galaxy clusters (DiMatteo. et al. 2004). Fortunately, all 
of these foregrounds have smooth continuum spectra whose 
contribution varies slowly with frequency. The contribution 
from these contaminants to the visibilities at the same base- 
line but two frequencies differing by Av will remain corre- 
lated for large values of Av, whereas the contribution from 
the HI signal will be uncorrelated beyond something like 
1 MHz or even less depending on the value of U . It is, in 
principle, straightforward to fit the visibility correlation at 
large Av and remove any slowly varying component thereby 
separating the contaminants from the HI signal. 

Finally, we briefly discuss the level of visibility corre- 
lation detectable in a radio-interferometric observation. We 
consider an array of N antennas, the observations lasting a 
time duration t, with frequency channels of width Sv span- 
ning a total bandwidth B. It should be noted that the effect 
of a finite channel width 8v has not been included in our 
calculation which assumes infinite frequency resolution. This 
effect can be easily included by convolving our results for the 
visibility correlation with the frequency response function 
of a single channel. Preferably, 8v should be much smaller 
than the frequency separation at which the visibility correla- 
tion become uncorrelated. We use S to denote the frequency 
separation within which the visibilities are correlated, and 
beyond which they become uncorrelated. The rms. noise in 
the a single visibility correlation is (Thompson, Moran & 
Swenson 1986) 



V((vv*y) 



2k B T S y 



A 



ef 



1 

~6vt 



(24) 



where Tsys is the system temperature and A e f is the effec- 
tive area of a single antenna. The noise contribution will be 
reduced by a factor 1/y/No if we combine No independent 
samples of the visibility correlation. A possible observational 
strategy for a preliminary detection of the HI signal would 
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be to combine the visibility correlations at all baselines and 
frequency separations where there is a reasonable amount of 
signal. This gives iV D = [N(N -l)/2] (B/Sv) (S/Su) whereby 



V((vv*y) = 



2k B T S YsY 



A 



ef 



N(N-1)BS 



(25) 



Using values for the GMRT (at ~ 325 MHz), 
(2k B T S Ys/A e f) = 144 Jy, B = 16 MHz, there are 14 an- 
tennas within U < 1000 where the signal is strong and 
S = 0.5 MHz beyond which the signal is uncorre lated we fi nd 
that it is possible to achieve noise levels of \/((VV*) 2 } — 
1.0 x 10 -9 Jy 2 which is below the signal with 200 hrs of in- 
tegration. 

We shall investigate in detail issues related to detecting 
the HI signal, the necessary integration time for different 
telescopes and signal extraction in a forthcoming paper. 
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3njj / h p c 2 Aioa 2 . dr(v) 
dv 
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where rim is the number density of HI atoms. In presence 
of peculiar velocity r(y) is given by 

K") = T ^n-, (A6) 



a?H{a) 



!/ e (l-«/c) 

where v is the line of sight component of the peculiar velocity 
of the HI. This gives us 

, dr(y) . A e r , 1 dv 



dv 



a 2 H(z 



rl 1 



aH{z) dr v 



(A7) 



where r„ is the same as r(y) without the effect of peculiar 
velocities, and where we have retained only the most dom- 
inant term in the peculiar velocity. Using this it is possible 
to write the optical depth as 



4.0 mK 

T 3 

PHI 
PH 



1 



0.02 

(1 + z) dv 



0/7 
h 



H 
H(z) 



(1 



H(z) dr v 



(A8) 



where pm/pH is the ratio of the neutral hydrogen to the 
mean hydrogen density. 



APPENDIX A: THE OPTICAL DEPTH OF THE 
REDSHIFTED HI 21 CM LINE. 



The optical depth is defined as 

r„(0= / a„{l')dl' 
Jio 



(Al) 



where a v is the absorption coefficient for the redshifted 21 
cm line given by 

hpfe 



An 



(u)[n Boi — niBio] 



(A2) 



Here <j>{v) is the line profile function, and Boi and Bio arc 
the Einstein coefficients (Rybicki. & Lightman 1979) for the 
21 cm line. We use r{y) to denote the comoving distance 
from which the 1420 MHz radiation will be redshifted to the 
frequency v for an observer at present. The HI radiation 
from a physical line element Al at the epoch when the HI 
radiation originated will be spread over a frequency range 
Av = (dr(v)fdv)^ 1 (Al/a) at the observer. Here a — 1/(1 + 
z) is the scale factor. The line profile function <f>(v), defined 
at the epoch when the HI radiation originated, is a step 
function of width aj Av. It then follows that 



(y) dl = a 



dr(v) 
dv 



and the optical depth can be written as 

-H)Boi(l- 

gi m 



4tt 



ffo no 2 , dr(v) 
av 



(A3) 



(A4) 



The Einstein coefficients are related by Boi = 3Bio = 
2jr~Aw , where Aio = 2.85 x 10~15 s _1 is the spontaneous 
emission coefficient of the 21 cm line . So the optical depth 
finally becomes 
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